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Abstract 

In this work, simple exact results are presented for summations 
in two-particle potential with long-range interactions. Polygamma 
function is used to evaluate summations. Results are found when a 
periodic media is consider. Periodic boundary conditions are applied 
by symmetric repetitions of a central cell. Contributions over all space 
are used in the problem. The potential depends strongly on the size 
of the system; however, forces are convergents anywhere. 

Introduction 

The thermodynamical extensivity imposes short-range interactions in classi- 
cal systems. Standard theoretical and experimental behaviors have been dis- 
cussed in details since several decades [TJ E] . However, anomalous behavior 
is obtained when the potential attractive tail behaves very slow. In recent 
years, much attention has been paid to physical systems with microscopic 
long-range interactions (i.e., see Hj and references therein). 

The most frequently applied way to discuss this type of systems is the 
Ewald method jHj, where neutralizing counter charges are introduced to en- 
sure convergence of the energy between one particle and the infinite replica- 
tions of other particle. A more recent approach is the Lekner method jSllZl, 
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where the symmetry of the lattice of identical computational cells is used to 
guarantee the convergence of the relevant interaction force. Here, no neutral- 
izing charges are considered and therefore, we adopt a Lekner-like procedure 
to evaluate the resulting interactions. We get in all cases (short and long- 
range interactions) convergent thermodynamical quantities. We discuss the 
most appropriate size of the system to evaluate such quantity. 

We have a central computational cell in one dimension with normalized 
size L — 1, this is — 1/2 < z < 1/2 where z is the variable of the position. 
This work is an analytical extension of a computational one introduced in [S] . 
Systems in one dimension are very important due to effects about ferromag- 
netism in one-dimensional monoatomic metal chains which has been reported 
recently [H]. Periodic boundary conditions are computed by repetition of the 
central cell to infinity. Particles interact with a potential given by 

oo 

v(z) = -J E g(\k + z\), (1) 

k=— oo 

where Jo is the hopping parameter and the summation on k represents all 
contributions over replicated images. Certainly, forces are obtained by f(z) = 
—zdv(z)/dz. 

Now, summations can be written in the following manner 

oo oo oo 

E g(\k + z\) = Y J g{k-z) + g {\z\) + Y J g{k + z). (2) 

k=— oo fc=l fc=l 

In the present work, the following expressions for potentials are taken into 
account 

, v _ J l/\z\ a power law interaction , , 

\ log(|^|) a logarithmic interaction. 

Replacing Eq.Q in Eq.(JTJ, v(z) is shown to be divergent; then, it is necessary 
to propose an alternative way for discussing this kind of potential. In present 
calculations, a particular summation (of l/{k+z) on k = 1, 2, 3, ...) is involved 
in the computation of potentials or forces . Let us use the definition of the 
digamma function [TU] 

where 7 = — .F(O) = 0.577215664901... that has been computed by several 
methods and accuracy IT2*] . We note that summations are divergents 
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what makes necessary to take a cutoff Now, we are going to approximate the 
summation from Eq.(J2J) by using the following approximated definition 



M I 

E — z = hm) - m (5) 

where M — > oo. We verify from Eq.Q) that if z = M is an integer, we can 
get 

M i 

Hm) = Y.t-i (6) 

fe=l K 

which converges only for finite values of M and diverges when M is infinite. 
In this way, M provides a natural cutoff for divergent summations. 



1 1/z potential 

From Eq.(^), the potential can be written as follows 

/ 1 00 1 00 1 \ 

and we find a close expression for the two-particle potential with a cutoff 
given by M replications in each direction of the lattice. If we replace Eq.© 
in Eq.(|7J), we can obtain an exact expression for the potential 

v(z) = -J Q + 2JF(M) - T{z) - F(-zj) . (8) 

from this equation, we can see that the potential v(z) depends on M. To 
eliminate that dependency, we will proceed to scale the potential as it has 
been conjectured [13]. For example, in this particular case, we can define a 
new scaled variable as follows 

/ x v(z) 
u(z) = — ^ 9 

This kind of scaling preserves the Legendre structure and it has been largely 
discussed in recent literature (i.e., see jSJCUlIS] and references therein). In 
Figure 1(a) we depict the v(z)/ Jq as a function of z. In this way, from Eq.© 
we verified that the potential does not saturate when M increases. This is a 
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typical behavior of a system with long-range interactions. Now, in figure 1(b) 
is depicted the scaled potential u(z)/J as a function of z, and we observe 
if M increases, the variable converges. The curve for M=l (dashed line) is 
far from the curve for M — 10 (long-dashed line), but as it is expected the 
curve for M = 10 is close to M — 100 (solid line), and we note a tendency 
to saturate when M increases. 

In addition to this, we can obtain the force between two particles and its 
replicated images as follows 

/(*) = z-^v(z) = zJ (1 + T\z) - , (10) 

where 

T { - m \z) = d m T{z)ldz m . (11) 

Polygamma function is obtained by differentiating the digamma function 
repeatedly. In the force calculation from Eq. fjlOj) . no dependency on M is 
obtained; this was appointed in a previous computational test [8 . In general, 
we use the symmetry of the lattice and contributions on a particle from 
images of the other one. 

The simplest case is one particle in the central cell and the potential must 
consider the potential on the central particle from its replicated images. This 
is therefore explicitly given by the following sum, 

m 

v = ~ 2J oJ2^ (12) 
k=i K 

and using Eq.© in the potential v, we can get 

v = -2J (F(M)+ 1 ). (13) 

Also, it is possible to see that the force is zero, this is / = kdv/dk = 0, for 
the simple case of one central particle in a linear system that only interacts 
with its images. 

2 Logarithmic potential 

When the potential is logarithmic, the force is wrote as follows 

oo-r / -i oo 1 oo-r\ 

/l * , -**£.^- , *t + £s+i-£i^ (14) 
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By using the identity 



1 00 1 



we can obtain an additional relation for representing an exact form for the 
force when the present logarithmic potential is considered, this is 

f(z) = zJ n cot(nz) = zJ a Q - T{z) + T(-z)\ . (16) 

The advantage for the second one is the explicit contribution to the force. 
The first term is the interaction on one particle due to the other particle in 
a particular cell, the second (and third) term is the contribution due to the 
images of the second particle situated on the right (left) of the central cell. 



3 Discussion 

On one hand, we have obtained a close expression for summations with the 
long-range interactions. Polygamma function is used and evaluation is fast 
and rapidly convergent. Potentials are strongly dependent on the size of 
systems. From Eq. fjlOj) . it is shown that the forces do not depend on the size 
of the system, such as it was obtained in a previous numerical work[8j. A 
direct cutoff is provided for potential summations by number M of repetitions 
of the central cell over all space. 

On the other hand, in the short range interaction regime when the at- 
tractive tail falls very fast as l/z 1+m law, the potential is given by 

00 I 

E n m — ( 17 ) 

According to Eq. fjllj) we can obtain an explicit definition for polygamma 
function, given by 

00 1 

where z is a real variable. This relation allows us a general expression for 
potential with the present power law short-range interactions. 

(_-\\m+l f(_-\\m+l \ 

v(z) = -Jo [ -^- (^r- + ^ m \z) + , (19) 
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being m = 1,2,3,.... In this result is shown that v(z) does not depend 
on size of the system when interactions are short ranged. Different terms 
are in relation to several contribution due to inter-particle interaction and 
contributions due to replicated images in each directions. 

Finally, this kind of generalizations are added to several contribution in 
this area, and the problem related to periodic boundary conditions is a topic 
for several studies. Thermodynamical limit for infinite systems with finite 
range of interactions is very well defined. However, as it has been shown, if 
the range of interactions is greater than the size of the system, thermodinam- 
ical quantities depend strongly on the size of system. It is very important 
to understand the microscopic behavior of the particles because the average 
from computer simulations comes from the average of finite number of par- 
ticles in systems. Advances on this topic will give us a view to understand 
the thermodynamical behavior for this kind of systems. 
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FIGURE CAPTION 
Figure 1 (a) Potential does not converge when M increases, (b) Scaled 
potential shows a tendency to saturate when M increases. Curves in the 
figure are depicted for M — 1 (dashed line) , M — 10 (long-dashed line) and 
M = 100 (solid line). 
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